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The statistical properties of random lattice knots, the topology of which is 
determined by the algebraic topological Jones-Kauffman invariants has been 
studied by analytical and numerical methods. The Kauffman polynomial in- 
variant of a random knot diagram has been represented by a partition function 
of the Potts model with a random configuration of ferro- and antiferromag- 
netic bonds, which allowed the probability distribution of the random dense 
knots on a flat square lattice over topological classes to be studied. A topolog- 
ical class is characterized by the highest power of the Kauffman polynomial 
invariant and interpreted as the free energy of a g-component Potts spin sys- 
tem for q — ► 00. It is shown that the highest power of the Kauffman invariant 
is correlated with the minimum energy of the corresponding Potts spin sys- 
tem. The probability of the lattice knot distribution over topological classes 
was studied by the method of transfer matrices, depending on the type of 
local crossings and the size of the flat knot diagram. The obtained results are 
compared to the probability distribution of the minimum energy of a Potts 
system with random ferro- and antiferromagnetic bonds. 
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I. INTRODUCTION 

New interesting problems are formulated, as a rule, in the boundary regions between 
traditional fields. This is clearly illustrated by the statistical physics of macromolecules, 
which arose due to the interpenetration of the solid state physics, statistical physics, and 
biophysics. Another example of a new, currently forming field is offered by the statistical 
topology born due to merging of the statistical physics, theory of integrable systems, alge- 
braic topology, and group theory. The scope of the statistical topology includes, on the one 
hand, mathematical problems involved in the construction of topological invariants of knots 
and links based on some solvable models and, on the other hand, the physical problems 
related to determination of the entropy of random knots and links. In what follows, we pay 
attention predominantly to problems of the latter kind that can be rather conventionally 
separated into a subfield of "probabilistic-topological" problems [0. Let us dwell on this 
class of problems in more detail. 

Consider a lattice embedded into a three-dimensional space, for which Q is the ensemble 
of all possible closed self-nonintersecting trajectories of iV steps with a fixed point. Let 
u denote some particular realization of such a trajectory. The task is to calculate the 
probability P for a knot of specific topological class for all trajectories to G Q. This can be 
formally expressed as 



2 



P{Inv} 



1 



5{lnv(c<j) — Invj 



(1) 



where Inv(cj) is a functional representation of the invariant for a knot on trajectory u, Inv 
is a particular topological invariant characterizing the given topological type of the knot or 
link, J\f(Q) is the total number of trajectories, and 8 is the delta function. The probability 
under consideration should obey the usual normalization condition P{Inv} = 1. The 



Based on the above definitions, it is easy to note that the probabilistic-topological problems 
are similar to those encountered in the physics of disordered systems and sometimes, as 
demonstrated below, of the thermodynamics of spin glasses. Indeed, the topological state 
plays the role of a "quenched disorder" and the functional P{Inv} is averaged over the 
trajectory fluctuations at a fixed "quenched topological state" , which is similar to calculation 
of a partition function for a spin system with "quenched disorder" in the coupling constant. 
In the context of this analogy, a question arises as to whether the concepts and methods 
developed over many years of investigation into the class of disordered statistical systems 
can be transferred to the class of probabilistic-topological problems, in particular, whether 
the concept of self-averaging is applicable to the knot entropy S. 

The main difference between the systems with topological disorder and the standard spin 
systems with disorder in the coupling constant is a strongly nonlocal character in the first 
case: a topological state is determined only for the entire closed chain, is a "global" property 
of this chain, and is difficult to determine it for any arbitrary subsystem. Therefore, we may 
speak of the topology of a part of some closed chain only in a very rough approximation. 
Nevertheless, below we consider a lattice model featuring a unique relationship between the 
topological disorder and the disorder in the local coupling constant for a certain disordered 
spin system on the lattice. 

Every time when we deal with topological problems, there arises the task of classification 
of the topological states. A traditional topological invariant, known as the Gauss invariant, 
in inapplicable because this Abelian (commutative) characteristic takes into account only 
a cumulative effect of the entanglement, not reflecting the fact that the topological state 
depends significantly on the sequence in which a given entanglement was formed. For ex- 
ample, when some trial trajectory entangles with two (or more) obstacles, there may appear 
configurations linked with several obstacles simultaneously, while being not linked to any 
one of these obstacles separately. In this context, it is clear that most evident questions con- 
cerning calculation of the probability of a knot formation as a result of the random closing 
of the ends of a given trajectory cannot be solved in terms of the Gauss invariant because 
this characteristic is incomplete. 

A very useful method of the classification of knots was offered by a polynomial invariant 
introduced by Alexander in 1928. A breakthrough in this field took place in 1975-1976, when 
it was suggested to use the Alexander invariants for classification of the topological state 
of a closed random trajectory computer-simulated by the Monte Carlo method. The re- 
sults of these investigations showed that the Alexander polynomials, being a much stronger 
invariants compared to the Gauss integral, offer a convenient approach to the numerical 



{Inv} 



entropy S of a knot of the given topological type is defined as 



S{Inv} = lnP{Inv} 



(2) 
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investigation of the thermodynamic characteristics of random walks with topological con- 
straints. The statistical-topological approach developed in || has proved to be very fruitful: 
the main results gained by now were obtained using this method with subsequent modifica- 
tions. 

An alternative polynomial invariant for knots and links was suggested by Jones in [|3],[4] . 
This invariant was defined based on the investigation of the topological properties of braids 
Jones succeeded in finding a profound connection between the braid group relationships 
and the Yang-Baxter equations representing a necessary condition for commutativity of the 
transfer matrix J/J. It should be noted that neither the Alexander, Jones, and HOMFLY 
invariants, nor their various generalizations are complete. Nevertheless, these invariants 
are successfully employed in solving statistical problems. A clear geometric meaning of the 
Alexander and Jones invariants was provided by the results of Kauffman, who demonstrated 
that the Jones invariants are related to a partition function of the Potts spin model |J. Later 
Kauffman and Saleur showed that the Alexander invariants can be represented by a partition 



function in the free fermion model 10 



We employed the analogy between the Jones-Kauffman polynomial invariant and the 



partition function in the Potts model with ferro (f-) and antiferromagnetic (a-) links pT|JT2 



to study the statistical properties of knots. In particular, the method of transfer matrices 
developed in [13-15] was used to determine the probability P{fx} of finding a randomly 
generated knot K in a particular topological state characterized by the invariant fx- 

The idea of estimating P{fx} is briefly as follows: (i) the Jones-Kauffman invariant is 
represented as a partition function of the Potts model with disorder in the coupling constant; 
(ii) the thermodynamic characteristic of the ensemble of knots are calculated using the 
method of transfer matrices for the Potts spin system. 

The paper is arranged as follows. Section 2 describes the model, defines the Jones- 
Kauffman polynomial invariant, and shows a connection of this invariant to the Potts model. 
Section 3 presents the numerical methods employed and introduces the necessary supple- 
mentary constructions. The results of numerical modeling are discussed and the conclusions 
are formulated in Section 4. 



II. KNOTS ON LATTICES: MODEL AND DEFINITIONS 

Consider an ensemble of randomly generated dense knots on a lattice in a three- 
dimensional space. A knot is called "dense" if the string forming this knot is tightly fit 
to the lattice, not fluctuating in the space. In this case, the knots of various topological 
types possess no configuration entropy and the probability of formation of a knot belonging 
to a given topological type is determined only by a local topology of the system. 

Not taking into account the fluctuational degrees of freedom of the lattice sites, this 
model is obviously very simplified. Nevertheless, this approach adequately reflects a physical 
situation, for example, in the statistical physics of condensed (globular) state of polymer 
macromolecules. Virtually all presently existing models taking into account topological 
restrictions with respect to the admissible spatial configurations of polymer chains assume 
a small density of the polymer, that is, describe a situation far from the compact state of 
globules. A hypothesis of the essentially new globular phase of a ring knotless macromolecule 
(crumpled globule) existing at a large polymer density (i.e., in the state of dense knots) was 
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formulated in fl6f . Despite an indirect experimental evidence for validity of the crumpled 
globule hypothesis ||17|| , direct observation of such objects in real experiments or computer 
simulations encounters considerable technical difficulties. In this context, investigation of 
the distribution of randomly generated dense knots of various topological types may provide 
for a better understanding of the structure of the phase space of knotted polymers in the 
globular state. 

Analysis of a string configuration in the three-dimensional space is not very convenient 
for determining the topological type of a knot. The standard method consists in project- 
ing the knot onto a plane in the general position (with no more than two knot segments 
intersecting at each point of the plane) and determining which segment line passes above 
(over crossing) or below (undercrossing) in accordance with the knot topology. Such a pro- 
jection, featuring over- and under-crossing, is called the diagram of a knot. In what follows 
we deal only with statistics of the knot diagrams. Evidently, this description implies an ad- 
ditional simplification of the model, but we believe that (in the phase of dense knots under 
consideration) the additional restrictions are not very significant since, as noted above, the 
knot entropy contains no contribution due to the string fluctuations. 

Thus, let us consider a square lattice of the size N = L x L on a plane, which is rotated 
for convenience by the 7r/4 angle. The lattice is filled by a (dense) trajectory featuring 
intersections at each lattice site with clear indications of the threads passing above and 
below (see an example of the 3x3 lattice in Figfja). As can be readily verified, a dense 
trajectory on a lattice with an odd L is unique, that is, represents a knot with the topology 
unambiguously determined by the pattern of over- and undercrossing and the boundary 
conditions. In this case, the probability of realization of the knot belonging to a given 
topological type is determined by the distribution of over- and undercrossing in all lattice 
sites. The problem considered in this study is to describe the distribution of dense knots 
over topological classes for various kinds of the over- and undercrossing distributions. 

A. Reidemeister moves and definition of the kauffman invariant 

In solving any topological problem, the main point consists in comparing the knots. A 
knot diagram on the plane obey the following Reidemeister theorem Jl8|: Two knots in the 
tree- dimensional space can be continuously transformed into each other if and only if the 
diagram of one knot can be transformed into that of the other knot by a sequence of local 
transformations (moves) of types I, II, and III (see Fig.g). 

As can be seen in Fig.[| the Reidemeister move I leads to the formation of a singularity 
on the plane during continuous fastening of the loop; this move is forbidden for smooth 
trajectories on the plane. Two knots are called regularly isotopic if their plane diagrams 
can be transformed into each other by means of the Reidemeister moves II and III. When 
the mutual transformation of the knot diagrams requires using the Reidemeister moves of 
all three types, the knots are called ambiently isotopic. 

Consider a two-dimensional knot diagram as a graph in which all intersection points 
(vertices) are characterized by the order of crossing (over- and under-crossing). Then each 
intersection point belongs to one of the two possible crossing modes. Let a kth point of the 
graph be characterized by the variable e& acquiring the values ±1, depending on the mode 
of crossing. 
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Let us define the algebraic Kauffman invariant as a sum over all possible variants of 
splitting the diagram at the vertices. According to this, each splitting is ascribed a certain 
statistical weight by the following rule: a vertex with e = +1 is given the weight A or B 
corresponding to the horizontal or vertical splittings, respectively; for a vertex with e = — 1, 
the weight B corresponds to the horizontal splitting and the weight A, to the vertical 
splitting. This definition can be illustrated by the following scheme: 



e= +1 




Thus, there are 2 various microstates for the diagram of a knot possessing N vertices. Each 
state uj of the knot diagram represents a set of nonintersecting and self-nonintersectiong 
cycles. The manifold of all microstates is denoted by {uj} (below, the braces {. . .} will 
denote summing over these states). Let S(u) be the number of cycles for the microstate uj. 
Consider a partition function 

f KR = ]T rfSM-i^M^-'M (4) 
M 

where the sum is taken over all 2^ possible splittings of the diagram; l(u) and iV — 1(uj) are 
the numbers of vertices with the weights A and B, respectively, for a given set of splittings 
of the microstate uj. Kauffman derived the following statement |J: a polynomial Jkr of the 
variables A, B, and d representing a partition function ([|) is a topological invariant of the 
regularly isotopic knots, if and only if the parameters A, B, and d obey the relationships: 
AB = 1 and ABd + A 2 + B 2 = 0. A proof of this statement in []9| is based on verification of 
the invariance of the partition function f^R with respect to the Reidemeister moves II and 
III. (The invariant for all tree Reidemeister moves is defined below.) The latter relationships 
pose the following restrictions on the parameters A, B, and d in Eq. ([|): 

d = -A 2 - A- 2 [b) 
which imply that invariant (H) is a polynomial of the single variable A. 



B. A partition function of the Potts model as a bichromatic polynomial 

Consider an arbitrary flat graph containing N vertices. Let ith vertex be characterized 
by a spin variable Oi (1 < Oi < q) and each edge of the graph, connect ing the ith and jth 
spins (1 < {i, j} < N), by the coupling constant Jij . The energy of the Potts model is 
defied as 0] 
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where the sum over {i,j} is taken only for the adjacent spins connected by edges of the 
graph. Then, the partition function can be expressed as 

M \{iJ} 

where {} denotes summation over all possible spin states, and the sum over {i, j} is taken 
as indicated above. The last expression can be written in the following form: 



Z = J2 II (! + VijSfa, crj)) , v id = exp ) - 1 (6) 

W{m'} v 



A pair of adjacent spins i and j introduces into the product term a contribution equal to 
exp(^-) for <jj = Oj and a unity contribution for Oi ^ <jj. Let us perform a procedure 
according to the following rules to a given spin configuration and the corresponding graph: 

• an edge is removed from the graph if a contribution to the above product from the 
spins connected by this edge is unity; 

• an edge is retained in the graph if the contribution from the spins Oi and cij connected 
by this edge is exp (4^-) . This procedure ensures a one-to-one mapping between the 
spin configuration corresponding to a product term in the sum and the related set 
of graph components. 

Consider a graph G containing M edges and C connected components (an isolated vertex 
is considered as a separate component). Upon summing over all possible spin configurations 
and the corresponding subdivisions of the graph G, we may present the sum (|6]) in the 
following form: 

M 

Z = Y,V C Ii %i (7) 

{g} m 

M 

where {G} denotes summing over all graphs and n is the product over all edges of the 

{m} 

graph G. Note that expression (|7|) can be considered as an analytic continuation of the 
Potts spin system to non-integer and even complex q values. 

For Vij = v, expression (|7|) coincides with a well-known representation of the partition 
function of the Potts model in the form of a bichromatic polynomial (see, e.g., [0,0)). The 
same expression is involved in the correspondence between the Potts model and the model of 
correlated percolation over sites and bonds suggested by Fortuin and Kastelleyn ||20||, which 



serves a base for the Monte Carlo cluster algorithms developed by Swedsen and Wang [21 



and Wolff [221 for the Potts model. 



C. Kauffman invariant represented as a partition function for the Potts model 

The Kauffman invariant of a given knot can be represented as a partition function for 
the Potts model on a graph corresponding to an arbitrary plane diagram of this knot, but 
we restrict the consideration below to an analysis of the knots on lattices. 
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Let us rewrite the Kauffman invariant in the form of a partition function for the Potts 
model determined in the preceding section, with M denoting the lattice knot diagram (see 
Figljja). The auxiliary variables s k = ±1 describe the mode of the knot splitting in each kth 
lattice site, irrespective of the values of variables e k = ±1 in the same vertices: 



Sk = +1 ) ( s k 



Let uj = {si, s 2 , ■ ■ ■ , sn} be the set of all variables characterizing the mode of splitting 
in the lattice containing N intersections. The Kauffman invariant (|4]) can be written in the 
following form: 

f KR (e u e N ) = £(-A 2 - A- 2 ) 8 ^- 1 exp fin A £ e k s k ) (8) 
M V fe=i / 

Here, {uj} indicates summing over all values of the variables s k (i.e. over all modes of 
splitting thy lattice diagram M) and the variables e k characterize a particular realization of 
"quenched" disorder in the system. Now we can demonstrate that configurations obtained 
as a result of splitting the diagram are in a one-to-one correspondence to configurations of 
the Potts model on a dual lattice. 

Consider the Potts model lattice A corresponding to the lattice M (see Fig.[L|b) , where 
open circles indicate positions of the Potts spins). The b i; j edge of lattice A corresponds to 
the kth site of the lattice diagram M. Therefore, disorder in the vertices of diagram M set 
by the variables e k coincides with the disorder in edges b it j of the lattice A, that is, with 
the disorder in coupling constants. Let us define the disorder of fejj for the lattice A via the 
coupling constants at the corresponding kth site of the lattice diagram M: 

^ j — efc, if the bond is vertical 
1,3 1 e k , if the bond is horizontal 

It should be recalled that the definition of the Kauffman invariant @ is based on split- 
ting the lattice diagram M into polygons representing a system of closed densely packed 
nonintersecting contours (Fig.|I|b). For a given configuration of splitting the lattice diagram 
M and the corresponding dual lattice A, we take the following agreement: all edges of the 
lattice A not crossed by polygons of the lattice M are labeled. In Fig.[L|b, the labeled edges 
are indicated by dashed lines. All the remaining edges of A are unlabeled. In these terms, 
the exponent in the partition function in Eq. ([8]) can be written as follows: 

hor. vcr. hor. vcr. 

E s k e k = E s k e k + E s k e k = E s k e k + E s k e k + E s k e k + E s k e k = 

& mark nonmark mark mark nonmark nonmark ^10^ 

hor. vcr. hor. ver. V / 

E hi + E Kj ~ E hi - E hi = -E6« + 2 E hi 



where 



E hi + E hi = E b 



all 
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Let rriu be the number of labeled edges and C u be the number of connected components in 
the labeled graph uj with N p vertices (each vertex corresponding to a Potts spin). The Euler 
relationship for this graph is 

S(u) = 2C W + m u -N p 
Now we can readily transform Eq. (Rf) to 



N ,2C. 



f KR (A, {6y}) = {-A 2 - A- 2 )-^ J] £ {-A 2 - A- 2 )^ J[ (a 2 ^ {-A 2 - A" 2 )) 

all {^} mark 

(11) 

by using ([T0|) and the fact that N is odd. Comparing expressions flTT| ) and (|6]), we obtain 
the equality 

E (A* + A- 2 )* '" ft (A 2b - (~A 2 - A" 2 )) =£11 CL + a,)) (12) 

{G} mark CT 

in which the right-hand part coincides with a partition function of the Potts model repre- 
sented in the form of a bichromatic polynomial. From this we obtain 

v M = A 2bi 'i (-A 2 - A- 2 ) = -1 - A Abi 'j 



q=(A 2 + A' 2 ) 2 

Since the disorder constants may acquire only discrete values b id = ±1, we may write the 
following expression for the coupling constants J^f 

Jij = T In (l - (A 2 + A~ 2 ) A 2bi A = T In (-A 4 ^) (13) 

Thus, we arrive at the following statement. The topological Kauffman invariant Jkr{A) 
of the regularly isotopic knots on the lattice M can be represented in the form of a partition 
function for a two-dimensional Potts model on the corresponding dual lattice A: 

f KR (A, {b id }) = K (A, {b itj }) Z (q(A), {./,. ; (/,,,. A)}) (14) 

where 

K (A, {btj}) = (A 2 + A- 2 Y (Np+1) exp f - In A £ bj] (15) 

is a trivial factor independent of the Potts spins. Here, the partition function of the Potts 
model is 

Z(g(A),{J M (6 M ,A)}) = £exp J ^- A) Ma,.a l )] (16) 

M \{iJ} 1 J 

with the coupling constant Jij and the number of states q given by the formulas 

^f = ln(-A 4 ^), q=(A 2 + A~ 2 ) 2 (17) 
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and the variable bij expressing disorder on edges of the lattice A corresponding to the lattice 
M. A relationship between b^j and e& is defined in Eq. (|J). 

A specific feature of the partition function (O) is the existence of a relationship between 
the temperature T and the number of spin states q. For this reason, T and q cannot be 
considered as independent variables. Once a positive q value is fixed, the variable A can 



formally acquire the complex values according to Eqs. (|17D , where logarithm can be taken 
of a complex argument. The appearance of complex quantities in the partition function can 
be interpreted in two ways. On the one hand, this implies expansion of the domain of the 
partition function to the complex plane. On the other hand, the parameters T and J it j do 
not enter explicitly into an expression for the Kauffman invariant and, hence, their complex 
values do not require any special consideration. Below we will be interested mostly in the 
probability distribution of the maximum power in the polynomial invariant and, therefore, 
can digress from particular values of the variables A, T, and Jij. 

Thus, we have determined an invariant Jkr for the regularly isotopic knots, the diagrams 
of which are invariant with respect to Reidemeister moves II and II. In order to obtain an 
invariant J'kr for the ambiently isotopic knots with oriented diagram, the corresponding 
partition function has to be invariant with respect to the Reidemeister moves of all types. 
Let us characterize each oriented intersection by a variable c& = ±1 according to the following 
rule 



( a ) X °k 





In addition, we define the knot twisting Tw(uj) as a sum of the Ck = bk values taken over 
all intersections: Tw{uo) = J2°k The invariant fxi{oj) of the ambiently isotopic knots can 

k 

be expressed as follows ||: 

!kM = f KR (u;)(-A) 3Tw ^ (18) 

It should be noted that our boundary conditions imply that Ck = bk- As is known, the 
Kauffman polynomial invariant fKi{u;A) of the ambiently isotopic knots as a function of 
variable A is equivalent to the Jones polynomial invariant fj(u;x) of the variable x = A 4 . 
Now we can use formulas (|14]) and ([15]) to express the Jones invariant through a partition 
function of the Potts model. Let a partition function of the Potts model have the form 
Z(t;q) = H(E,q)t~ E where t = e E ; E denotes the energy levels over which the sum is 

E 

taken, and H(E, q) is the degree of degeneracy of an energy level E for a given q value. 
Taking into account that x = A 4 = — t and using formulas (0), (|15|) and (|T8|), we obtain 
the following expression for the Jones polynomial invariant: 

N p +l+ J2 h,j 

f J = J2 H ( E ^ = 2 + x + x~ l )(~x)- E (l + x)- N "- 1 (-^) <«> (19) 

E 

In what follows, by the maximum power of a polynomial invariant we imply the maximum 
power of variable x in the Jones polynomial invariant fj. It will be born in mind that the 
power of the Kauffman polynomial invariant fxi of the ambiently isotopic knots is obtained 
from the corresponding power of the Jones polynomial invariant through simply multiplying 
by a factor of four. 
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As indicated above, out task is to calculate the probability P{fj} of finding a knot on 
a lattice in the topological state with a preset J ones-Kauffman invariant fj(x, {e^}) among 
all the 2 N possible disorder realizations {e^}, k = 1, . . . , N This probability can formally be 
written as 



Thus, the topological disorder determined by a random independent selection of inter- 
sections of the e = +1 and e = — 1 types represents a random quenched external field. When 
we use the Jones invariant, each topological class (homotopic type) is characterized by a 
polynomial. In this case, precise identification of the homotopic type for a knot on a lattice 
containing N intersections will require N variables. Since the number of various homotopic 
types increases as 2^ , it is very difficult to study the probability of each separate homotopic 
type characterized by N variables. For this reason, we introduce a simplified characteristic 
of a knot on a lattice - the maximum power n of the polynomial invariant fj(x): 



For a trivial knot n = and, as the knot complexity grows, the maximum power increases 
(not exceeding N). Thus, ensemble of all of knots on a lattice can be divided into subclasses 
characterized by the power n (0 < n < N) of the Jones polynomial invariant. In these 
terms, we will study the probability that a randomly selected knot belongs to one of these 
subclasses (and is characterized by the maximum power n of the Jones invariant). Random 
knots can be generated by two methods: 

1. place a fixed number of intersections with e = — 1 on the knot diagram with N vertices 
(accordingly, the other vertices belong to intersections of the e = +1 type); 

2. place intersections of the e = — 1 and e& = — 1 types in each vertex with the probabil- 
ities p and 1 — p, respectively. 

As can be readily seen, a trivial knot having all intersections of the e = +1 type corresponds 
to a partition function with the ferro (f-) and antiferromagmetic (a-) bonds distributed in 
accordance with the rule @. An impurity (corresponding to = — 1 of the knot diagram) 
will be considered as a change in the sign of bij relative to the values characterizing trivial 
knots. Note that we must differentiate between the notion of impurity (a change in the sign 
of from the a-bond with b^j — — 1. For a trivial knot with all = 1, the lattice contains 
no impurities while containing the a-bonds. 

III. AUXILIARY CONSTRUCTIONS AND NUMERICAL METHODS 

A. The form of a lattice for the Potts model and the positions of ferro- and 

antiferromagnetic bonds 

Now we describe the geometry of a lattice for the Potts model corresponding to a knot 
diagram of the N = L x L size. This study is restricted to square lattices, although all 





(20) 
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considerations remain valid for rectangular lattices as well. Figure ||a shows an example of 

the trivial knot on a lattice with N = 5x5. Positions of the Potts spins corresponding to this 

lattice diagram are indicated by circles. The Potts spin lattice corresponding to this knot is 

depicted in Fig.||b, where the f-bonds (that would be horizontal in Fig.|l|) are indicated by 

solid lines and the a-bonds (that would be vertical in Fig.[l] ) are indicated by dashed lines. 

In FigFJc, the same lattice is transformed into a rectangular one of the Lh x L v type, where 

Lh = L + 1 and L v = (L + l)/2 (in this particular case, L v = 3, Lh = 6. The partition 

function of the Potts model is studied below for a rectangular lattice of this type. 

The fact that the height of the lattice for the Potts model is half that of the lattice 

knot diagram (see Figs.|3]a-c) is very convenient in the case of using the method of transfer 

matrices, where the computational time expenditure exponentially depends on the lattice 

height. Let N = L x L be the total number of bonds in the lattice, with iV + = $(bij, 1) 

{'•./} 

representing the number of f-bonds and iV_ = ${bij,— 1) the number of a-bonds. An 

{id} 

impurity, corresponding to = — 1 in the initial lattice, is described by a change in sign of 
the corresponding constant 6y for the Potts model lattice. 

It should be emphasized once again that by impurity we imply a change in the type of 
intersection from = +1 to = — 1, which corresponds to a change in sign of the coupling 
constant bij for this intersection, rather than to the a-bond (such bonds are present in the 
Potts model lattice of the trivial knot). A trivial knot is characterized by the absence of 
impurities on the Potts model lattice N + = AL = Note that the Kauffman in- 

variants for a knot and its mirror image are different. This is due to the following property 
of the Jones invariant: the invariant of the mirror image of a knot is obtained by substi- 
tuting t — > t^ 1 (and A — ► A^ 1 ), after which the distribution of powers in the polynomial is 
asymmetric with respect to the substitution p — > 1 — p. 



B. The method of transfer matrix 

It should be recalled that we are interested in determining the probability distribution 
of the maximum power n (averaged over various types of intersections in the knot diagram) 
of the Jones-Kauffman invariant. According to definition fl20|) , the coefficient at the term of 
the maximum power n in the Jones polynomial of a randomly generated knot is insignificant. 
Let us fix a certain value no, generate an ensemble of random knot diagrams, and determine 
the fraction of knots the polynomials of which contain the maximum power n = n . 

The traditional approach to numerical analysis of a g-component Potts spin system 
assumes that each possible state of the column L v spins corresponds to an eigenvector of 
the g L "-dimensional transfer matrix. With computers of the P-II-300 type and a reasonable 
computational time, the maximum possible size of the transfer matrix is 100-200. This 
corresponds to a spin strip width L v = 8 for q = 2 and L v = 4 for q = 3. Obviously, this 
method cannot be used for investigation of the Potts model with large q values. 

A method of transfer matrices applicable to the case of arbitrary q was developed by Blote 



and Nightingale |i3fl . The basic idea of this method is that each eigenvector corresponds 
to a subdivision of the column into clusters of "coupled" spins. By definition, the spins 
in each cluster are parallel. Therefore, the number of "colors" q can be considered simply 
as a parameter taking any values (including non-integer and even complex). A detailed 
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description of this method can be also found in [p4| , [15| , where this approach was used to 
study the Potts model with impurities on the bonds and to determine the roots of a partition 
function for the Potts model in the complex plane. 

Below we will briefly describe the principles of determining a basis set in the space of 
coupled spins and constructing transfer matrices in this basis. This description follows the 
results reported in flSfl . 

A basis set corresponds to a set of various subdivisions of the spin column into clusters 
of coupled spins. We must take into account that the subdivision should admit realization 
in the form of a flat graph on a half-plane. Now we can describe a recursive procedure for 
determining the basis set. Let a basis set be available for a column with the height L v . 
Upon adding one more spin from below, we obtain a column with the height L v . Then we 
take, for example, a subdivision V\(L V ) corresponding to the first basis vector L v and begin 
to generate subdivisions corresponding to the basis vectors for the column L v by attaching 
the added spin to those existing in v 2 (L v ). In each step, we check for the possibility of 
attaching this spin within the framework of the flat graph on a half-plane, after the attempt 
at attaching the added spin to all existing clusters, we generate a subdivision in which the 
added cluster is separate. Then we take the next subdivision (for example, V2(L V )) for the 
L v column and repeat the generation procedure. 

The procedure begins with a single spin as depicted in Fig.[|. This spin is assigned the 
index of unity as belonging to the first cluster. Then another spin is added from below to 
the same column. Accordingly, the strip of two spins may occur in one of the two states: 
v±(L v = 2), in which the spins are coupled and belong to cluster 1, or V2(L V = 2), in which 
the spins are not coupled and belong to clusters 1 and 2 (Fig||). A basis set for the strip 
of three spins is obtained by adding another spin from below and coupling this spin to 
all clusters. The clusters are enumerated by integers top to bottom, while vectors in the 
basis set are enumerated in the order of generation. As can be readily seen, the first basis 
vector corresponds to a state in which all spins are coupled to each other and belong to 
the same cluster, while the last basis vector corresponds to the state where all spins are 
uncoupled and the number of clusters equals to the number of spins in the column. The 
total number of vectors in the basis set for a strip of m spins is determined as the Catalan 
number C m = ( 2 ™J with the corresponding index number (see |15[| ). 

Let us denote by Vi(\k) a subdivision obtained upon separating the kth spin from the 
subdivision V{ to form a separate cluster; Vi({k, /}) will denote a subdivision in which the 
clusters containing kth and Ith spins are combined to form a common cluster. Following a 
method described in ITSfl, let us form the matrix 



A j (*0 = <J( u i(\*0> u *) + q5(v j (\k),v j ) 

As can be seen, an element of the matrix D it j(k) is non-zero if the kth spin is not coupled 
to any other spin in the subdivision i>j. When the same spin is not coupled to any other in 
the jth subdivision as well, the matrix element is q, otherwise it is taken equal to unity. We 
will also introduce the matrix C(k,l) defined as 

C iJ (k,l) = S(v j ({k,l}),v i ) 

the elements of which are equal to unity if and only if the subdivision Vi is obtained from 
the subdivision Vj by combining clusters containing the kth and Ith spins. 



13 



Now let us consider a lattice structure with the first two rows such as depicted in Fig.[|d 
(for an even index j + 1 in the second column) and Fig.|5|e (for an odd index j '• + 1 in the 
second column). The spin columns are enumerated by index j: 1 < j < L h , while the spin 
position in a column is determined by index i: 1 < % < L v . 

According to these constructions, the transfer matrices of the type T even (j + 1 = 2k), 
which corresponds to adding an even column j + 1 = 2k after an odd one j — 2k — 1 (we 
add the columns from the right-hand side), contain bonds between the spins (Tij=2k-i an d 
°~i+i,j=2k- The transfer matrices of the type T odd (2j + 1), which corresponds to adding an 
odd column j + 1 = 2k + 1 after an even one j = 2k, contain bonds between the spins o"j J=2 / c 
and (Ti-i J j+i=2k+i- 

Let the contribution to a statistical weight corresponding to the f-bond b{pij , ay j/) = 1 
be expressed as 

t = exp(/3J) 

where (3 = ~. Then the a-impurity b(cFij , ji) = — 1 corresponds to the weight 

r 1 = exp(-0J) 

and the total contribution of a given configuration to the statistical weight of the system is 
(WWj'), Let us define a matrix for the horizontal bond ((Tij, a s 

P(i,j) = I(t b ^^-l) + D(i,j) 

and the matrices for "sloped" bonds as 

R even (iJ) = I + C{i,]){t b ^^+^+^ - 1) 
R° dd (i,j) = I + (7(i, j)(t 6 K;,^-u+i) _ i) 

where I is the unit matrix. 

Then the transfer matrices can be expressed as follows: 

T even^ = P (L V , j)R^ n (L v - I, j)P(L v - 1, j) . . . R even (l, j)P(l, j) 

T° dd (j) = P(l,j)R° dd (2,j)P(2,j)...R° dd (L v ,j)P(L v ,j) 
and the partition function is 

W2-1 

Z(t)=u T T™ en (L h )(T odd (2k + l)T even (2k))v max 

fc=i 

where v max is a basis vector with the maximum number corresponding in our rep- 
resentation to the state in which all spins belong to different clusters and u T = 
{ q NcL(vi)^ q NcL{v2)^ ^ q N C L(v m ax)} go t h at yT v . = q N CL (vi) ? ^qj-q N CL {vi) v& the number of 

various clusters in the subdivision V{. For the first basis vector, Ncl{v\) = 1 (a single 
cluster) and the last vector corresponds to N C Liy max ) = L v - L v different clusters. 
The general algorithm of the calculation is as follows. 

1. Generate basis set vectors for a spin strip of the required width; 

2. Use this basis set to generate the matrices R(i,j), P(i,j) in cases of f- and a-bonds 
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3. Generate a distribution of impurities on the lattice using a random number generator; 

4. Generate the transfer matrices for the Potts spin model and the polynomial invariant; 

5. Calculate a partition function for the Potts model, the polynomial invariant, the min- 
imum energy, and the maximum power of the polynomial. 

6. Repeat points 3-5 and perform averaging over various realizations of the distribution 
of impurities on the lattice. 

The error was determined upon averaging over ten various series of calculations. For 
each realization of the impurity distribution, the program simultaneously determines both 
the polynomial invariant and the partition function of the Potts model for an arbitrary (but 
fixed) q value. This method allows us to study correlations between the maximum power of 
the polynomial and the minimum energy within the framework of the Potts model on the 
corresponding lattice. 

IV. RESULTS OF CALCULATIONS 

A. Correlations between the maximum power of Jones polynomial of the lattice knot 

and the minimum energy of the Potts model 

Let us consider dependence of the maximum power of the Jones polynomial on the type 
of the partition function for the Potts model. Formula flT9| ) explicitly relating the Jones- 
Kauffman invariant to the Potts model shows that, if the variable x were not entering into an 
expression for the degeneracy of the energy level H(E; q = 2 + x + x -1 ), the maximum power 
of the polynomial invariant would always correspond to a term of the partition function with 
the minimum energy. 

Taking into account dependence of the degree of degeneracy on the variable x, we can see 
that contributions to the coefficient at this power in some cases mutually cancel each other 
and this coefficient turns zero. A simple example is offered by a system free of impurities, 
in which case the minimum energy is E min = —(N + l)/2 and the Jones polynomial is 
identically unity (with the maximum power being zero). 

Nevertheless, there is a strong correlation between the maximum power n of the Jones 
polynomial and the minimum energy E m i n of the corresponding Potts model. Since the 
minimum energy E min of the Potts model is always sign-definite and cannot be positive, we 
use below a positive quantity representing the absolute value of the minimum energy \E min \ = 
—E min . Figure ^|a shows a joint probability distribution P (j^, ^'"^ of the normalized 

maximum power ^ and the normalized minimum energy ^pJ. obtained for a lattice with 
N = 49 and an impurity concentration of p = 0.5 by averaging over Nl = 10 5 impurity 
configurations. Here and below we use only the normalized quantities determined on the 
[0, 1] segment for the energy and [—1,1] for the maximum power of the polynomial (the 
latter value can be negative). This allows us to plot the curves for various lattice dimensions 
on the same figure. In Fig.|B|a, the probability distribution is described by the level curves, 
the spacing between which corresponds to a probability difference of 0.001. As can be seen, 
there is a strong correlation between the maximum power of the polynomial invariant and the 
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minimum energy of the corresponding Potts model. This relationship can be quantitatively 
characterized by the the coefficient of correlation. 

It should be recalled that the coefficient of correlation between random quantities x\ and 
X\ with mathematical expectations (x\) and (x 2 ) and dispersions Ax\ = (xf) — {x\) 2 and 
A^2 = (x 2 ,) — (X2) 2 , (where (. . .) denotes averaging) is determined by the formula 



corr(xi, x 2 ) 



The correlation coefficient equal to ±1 corresponds to a linear relationship between x% and 
x 2 . 

The values of the coefficient of correlation between the maximum power n max of the 
polynomial invariant and the minimum energy |-£7 m j n | for the impurity concentration p = 0.5 
on lattices with various the linear dimensions L = 3,5, 7, 9, and 11 are presented in Table 1. 
These values were obtained by averaging over Nl = 10 5 impurity realizations for L = 3 4- 7, 
Nl = 2 x 10 4 for L = 9, and Nl = 10 3 for L — 11. As is seen, the correlation increases with 
the lattice size. 



TABLE 1 . Mean value of the correlation coefficient corr (n 
statistical error for several lattice sizes. 



max ; l-^min I 



and the corresponding 



L 


corr (n, \E min \) 


Acorr (n, \ E min \) 


3 


0.4871 


0.0021 


5 


0.6435 


0.0022 


7 


0.7205 


0.0007 


9 


0.7692 


0.0013 


11 


0.7767 


0.0129 



n \E„ 
.A" 



Figure |6|b shows the approximation of these data by a power function corr (- N , v 

1.04(4) — 1.126(14)L _0 - 65< ^ 7 \ Naturally, the correlation coefficient cannot be greater than 
unity. The last expression is the result of ignoring the higher terms of expansion in powers 
of jk The approximation shows that there is a relationship between the maximum power 
of the polynomial invariant and the minimum energy of the Potts model, the degree of 
correlation increasing with the lattice size. In what follows, the results for the maximum 
power of the polynomial invariant will be accompanied by data for the minimum energy. 

The presence of a correlation between the maximum power n of the Kauffman polynomial 
invariant and the minimum energy E min of the corresponding Potts model for each particular 
impurity realization on the lattice is of interest from theoretical standpoint and can be used 
in the numerical experiments. At present, the limiting lattice size Nl used in the method 
of transfer matrix does not exceed Nl = 11, which is determined by a large volume of 
necessary computations. At the same time, the minimum energy for the Potts model can 
be calculated within the framework of the standard Monte-Carlo algorithm, which poses a 
much smaller requirements to the computational facilities and, hence, can be readily applied 
to the lattices of significantly greater size. 
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B. The probability distribution of the maximum power of the polynomial invariant 
and the minimum energy of the corresponding Potts model 

Here we present the results of determining the probability distribution P(n) of the max- 
imum polynomial power n. Figure |[]a shows the P(n) curves for L = 7 and the impurity 
concentrations p = 0.1 (depicted by crosses), 0.2 (squares), and 0.5 (circles). The data were 
obtained by averaging over Nl = 10 5 impurity realizations. The statistical errors are smaller 
than the size of symbols. 

As can be seen, the P {^\ curve for small impurity concentrations is nonmonotonic. 
As the p value increases, the probability distribution function becomes monotonic and ap- 
proaches in shape to the Gauss function. Figure shows the corresponding probability 
distribution of the normalized minimum energy ™' n obtained for the same lattice size 
(L = 7) and impurity concentrations (denoted by the same symbols) . This function appears 
as more monotonic and approaches the Gauss function already at small impurity concentra- 
tions (p ~ 0.2). 

The shapes of the probability distribution observed for a fixed impurity concentration 
(p = 0.5) and various lattice dimensions are shown in |TT| a (for the maximum polynomial 
power) and Fig.|TT|p (for the minimum energy). Data for the lattice size L = 3 (squares), 
^(squares), and 7 (circles), and 9 (triangles) were averaged over Nl = 10 5 (L = 3, 5, and 7) 
and data for L = 9, over Nl = 5 x 10 3 impurity realizations. The probability was normalized 
to unity: P{n) = 1 For this reason, an increase in the lattice size is accompanied by 

n 

growing number of the values that can be adopted by the normalized maximum power of 
P{jr)i while the value of the probability distribution decreases approximately as i. As 

is seen, the probability distribution P for lattices of smaller size (L = 3 and 5) is 
nonmonotonic; as the L value increases, the distribution becomes a smooth Gaussian-like 
function. 

Thus, we may conclude that the probability distribution P (^J of the maximum power 
of the polynomial invariant for small-size lattices is nonmonotonic as a result of the bound- 
ary effects even for a considerable impurity concentration. As the lattice size increases, 
the nonmonotonic character disappears and the probability distribution becomes a smooth 
function. Apparently, we may ascertain that, whatever small is the impurity concentration 
p, there is a lattice size N (N ^> i) such that the corresponding probability distribution 
is smooth and Gaussian-like. Thus, we may suggest that the probability distribution P(n) 
on the lattices of large size is determined by the Gauss function, the main parameters of 
which are the mathematical expectation (mean value) of the maximum polynomial power 
and the dispersion. Plots of the mean value of the maximum polynomial power and the 
corresponding dispersion W% nat as functions of the impurity concentration p are presented 
in Figs. ||a and Fig|| respectively, for the lattices size of L = 3 (crosses), 5 (squares), and 
7 (circles). The data for each point were obtained by averaging over Nl = 10 5 impurity 
realizations for L = 3 and 5 and Nl = 5 x 10 4 for L = 7. 

As expected, both the mean value and the dispersion of the maximum polynomial power 
turn zero for p = and 1 (trivial knot) and reach maximum at p ~ 0.5. Note that these 
functions are not symmetric with respect to the transformation p —>■ 1 — p, since the Jones 
polynomial invariant of a mirror knot (with all overcrossing changed for undercrossing and 
vice versa) is obtained by changing variables x — > x^ 1 , whereby the maximum polynomial 
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power of the mirror knot corresponds to the minimum power of the original polynomial, 
taken with the minus sign. However, this asymmetry disappears with increasing lattice size 
N as a result of increase in the amount of impurities and in the number of possible impurity 
realizations employed in the averaging. For the comparison, Fig.|9|b shows the function 
2^(1 — p) representing the dispersion of the distribution function in the hypothetical case 
when the maximum polynomial power n is a linear function of the number of impurities M. 
Thus, a difference between W% not and — p) characterizes dispersion of the distribution 
of the maximum polynomial power at a fixed number of impurities M (it should be recalled 
that the impurity occupies each lattice site with a probability p and the total number of 
impurities M fluctuates). 

Figures |lOl a and |lO| b (with the parameters and notations analogous to those in Fig.|a 
and Fig.^b) shows data for the absolute value of the mean minimum energy / \ an d 
the corresponding dispersion Wp otts . The asymmetry of the mean minimum energy plot is 
related to the fact that the number of f-bonds is greater than that of the a-bonds by one at 
p = and is smaller by one, at p — 1. The probability distribution of the minimum energy 
is treated in more detail in Subsection |A2| of the Appendix. We have studied dependence of 
the mean normalized maximum polynomial power on the lattice size L for p = 0.5 (Fig.|TT|a). 
The results were averaged over Nl = 10 5 impurity realizations for L = 3,5, and 7 and over 
Nl = 1.5 x 10 4 realizations for L = 9. As can be seen, the (jj) value tends to a certain 
limit with increasing L — > oo . We approximated the results by a power function to obtain 

~ 0.334(8) - 0.41(2)L- a48(5) (21) 

Analogous data for the normalized minimum energy in the Potts model are depicted in 
Fig.|i~l"|b. The mean absolute value of the minimum energy decreases with increasing lattice 



size and can be approximated by a power function of the type 

^Knr^ ^ _ .4 18 5( 7 ) + .119(3)L- L11(4 ) 

Some features of the minimum energy distribution in the Potts model with random ferro- 
and antiferromagnetic bonds are treated in the Appendix for q > 4, in which case analytical 
expressions can be obtained for small (p ~ 0) and large (p ~ 1) impurity concentrations. 



V. CONCLUSIONS 

The results of this investigation can be formulated as follows. 

1. An analysis of the relationship between the partition function of the Potts model 
and the Jones-Kauffman polynomial invariant, in combination with the results of numerical 
calculations for the Potts model, allowed us to study the probability distribution P(j^;p, N) 
of the maximum polynomial power n for various concentrations (probabilities) p of an "im- 
purity" representing intersections of the e = — 1 type in the sites of a square lattice with 
AT = L x L for L = 3 - 11. 

2. For the lattices of small size with small impurity concentrations p, the probability 
distribution P(^;p,N) of the normalized maximum polynomial invariant power n is not 
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smooth, showing an alternation of more and less probable states. This behavior indicates 
that the knots of certain topological types are difficult to realize for p < 1 on a square 
lattice result of geometric limitations. 

3. As the impurity concentration (probability) p increases, the probability distribution 
P(jr',p, N) of the normalized maximum polynomial invariant power n becomes smooth, with 
the shape approaching that of the Gauss function. A typical value of the "knot complexity" 
rj for p = 0.5 and sufficiently large lattices (N ^> 1) can be obtained by extrapolating the 
expression fl2~l~|): 



4. There is a correlation between the maximum power of the polynomial invariant of a 
knot and the minimum energy of the corresponding Potts lattice model for q > 4. 

5. An analytical expression was obtained for the probability distribution function 
P{ j$ n ',p, N) at relatively small p ~ i impurity concentrations. 

We have studied the knots with a square lattice diagram. However, all calculations can 
be generalized to the case of rectangular lattices. In connection with this, it would be of 
interest to check whether the distribution of knots over the topological types only depends 
on the number of intersections on the lattice diagram or it depends on the diagram shape 
as well. We may expect that the probability distribution function for a strongly elongated 
rectangular diagram would differ from the distribution for a square diagram with the same 
number of intersections. 

The method developed in this study is applicable, in principle, to the investigation of 
knots with arbitrary diagrams (including the case of knots with the diagrams not tightly fit 
to a rectangular lattice), provided that a configuration of the Potts system corresponding 
to this diagram is known. 

We believe that the proposed combination of analytical and numerical methods for the 
investigation of topological problems using the models of statistical physics offers both a 
promising means of solving such topological problems and a new approach to the standard 
methods of investigation of disordered systems. This can be illustrated by the following 
fact. One of the main concepts in the statistical physics is the principle of additivity of the 
free energy of a system, that is, proportionality of the free energy to the system volume 
N. By interpreting the free energy as a topological characteristic of the "complexity" of 
a knot, we may conclude that the complexity of a typical knot increases linearly with the 
system volume N. This property is, in turn, well known in the topology (outside the 
context of statistical physics), being a fundamental manifestation of the non-Abelian (non- 
commutative) character of the phase space of knots. 
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APPENDIX: MINIMUM ENERGY DISTRIBUTION IN THE POTTS MODEL 
WITH RANDOM FERRO- AND ANTIFERROMAGNETIC BONDS 



1. Independence of the minimal energy upon q for q > 4 

The minimum energy of a spin system is independent of q for q > 4. In this study of the 
topological invariants of random lattice knots, we are interested in determining behavior of 
the system for x — > oo(see Eq. (POD)- With an allowance for the relationship q = 2+x+x~ l , 
this implies q — > oo for the corresponding Potts model. Some features of the dependence 
of the free energy of the Potts spin system on the number of states q can be established 
based on simple considerations. In cases when the lattice contains no impurities and the 
a-bonds are arranged as depicted in Figs.[3]b and |]c, the minimum energy is independent of 
q for q > 2. Under these conditions, two spin values a = {1,2} are sufficient to "create" 
a configuration corresponding to the minimum energy. In the presence of impurities, the 
minimum energy depends on q for q > 2. This is illustrated in Fig.|5|a for a 5 x 5 lattice with 
a single impurity on the bond between a\,\ and a\^, spins. For the spin configuration energy 
to reach a minimum value of E m i n = — 12 for the given distribution of bonds, it is necessary 
that all spin variables in the lattice sites connected by f-bonds J = 1 (in Fig.[5], these sites are 
indicated by open and filled circles) would acquire the same values, while the spin variables 
in the sites cr^i and a^ 2 connected by a-bonds J = — 1 (in Fig.[5]a, these sites are indicated by 
filled and hatched circles) would acquire different values (e.g., a 1; i = 2, a X 2 = 3). However, 
it is impossible to assign the values of spin variables in a system with q = 2 so that the spins 
in clusters cr^i, a\ j2 (as well as in the adjacent cluster) indicated by unlike symbols (black 
versus open or hatched) were different. For this reason, a minimum energy of the spin state 
for q = 2 is E min = —11 (instead of —12). 

As was demonstrated, a minimum energy of the spin configuration in the Potts model 
corresponding to a given distribution of a-bonds may depend on the number of spin states 
q. This fact can be represented as follows: to reach the state with minimum energy, it is 
necessary that the spin variables in the lattice sites connected by f-bonds would acquire 
for the most part the same values (so as to form clusters), while the spin variables in the 
sites connected by a-bonds would be possibly different. Thus, a given distribution of the 
a-bonds corresponds to a subdivision of the lattice into independent clusters of spins. Spins 
belonging to the same cluster are connected predominantly by the f-bonds, while spins of 
different clusters are connected by a-bonds. 

The energy reaches minimum if the adjacent clusters in a given subdivision possess 
different values of the spin variable. We may bring each value of the spin variable into 
correspondence with a certain color. Then a minimum energy corresponds to the lattice 
subdivision into clusters painted so that all spins in one cluster are of the same color, whereas 
adjacent clusters have different colors. Figure |5]a shows an example of the configuration 
which cannot be painted in this way using two colors, while three colors allow reaching the 
goal. In mathematics, there is a theorem concerning the task of "painting maps", according 
to which any configuration on a surface possessing a a topology of a sphere can be painted 
using four (or more) colors. In other words, any subdivision of a lattice into clusters can 
be painted using four (or more) colors so that the adjacent clusters would possess different 
colors. If each color corresponds to a certain value of the spin variable q = {1,2,3,4, ...}, 
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we may assign the q values (for q > 4) so that spins in the adjacent clusters would possess 
different values (colors). There is no impurity configuration (and the corresponding lattice 
subdivision into clusters) such that four values of the spin variable would be insufficient to 
reach the state of minimum energy. 

The above considerations allow us to formulate the following statement: For an arbitrary 
configuration of a-impurities on a Potts model lattice, a minimum energy of the spin system 
is independent of the number of spin states q for q > 4. 

As is known, the Potts model exhibits a first-order phase transition at q > 4. Note also 
that the parameter in the Kauffman invariant becomes real just for q > 4. 

2. Distribution of the minimum energy at small (p ~ 0) and large (p ~ 1) impurity 

concentrations 

Consider the Potts model on an L v x lattice with a total number of lattice sites N and 
q > 4. Let M be a fixed number of impurities of the types. In the absence of impurities 
(M = 0), the number of f-bonds is N + = and the number of a-bonds is N_ = ^p-. 
The arrangement of a-bonds is illustrated in Fig.|5|a. The minimum energy in the absence of 
impurities is E min = —N + . 

TABLE 2. The probability distribution of energy P(E m i n ; M, N) for a fixed number of impu- 
rities M = 0, 1, 2>iV - 2, N - 1, N on a lattice with iVsites (C = \^N). 



I 


II 


III 


IV 


M 


Link type 


Emin 


P{E min , M) 







N + l 
2 


1 


1 


+ 


2 r x 


N+l 
2N 


1 




N + l 


N-l 
2N 


2 


++ 


N + x + 2 


(N+1)(N-1) 
4N(N-1) 


2 


+- 


N+l + 1 


(N+1)(N-1)-W(C-1) 


2 ~ ± 


2N(N-l) 


2 


+-&— 


N+l 


(N-l)(N-3)+W(C-l) 


2 


4NIN-1) 


2 




N+l i 
2 1 


4(0-2) 
N(N-l) 


2 




N+l 9 
2 Z 


4 

N(N-l) 


N-2 


++ 


-^i + 2 


(N-l)(N-3) 
4N(N-1) 


N-2 


+- 




(JV-l)(JV-3)-16(6'-l) 


-^ + 1 


2JV(JV-1) 


N-2 


+-&— 


N-l 


(JV+l)(JV-l)+16(C-l)-4 


2 


4N(N-1) 


N-2 


+-&— 


N-l i 


2(N-1)+2(C-1)) 


2 1 


N(N-l) 


N-2 




N-l o 

2 Z 


^ 


N-l 


+ 


N-l + 1 
2 r x 


2W 


N-l 




N-l 
2 


7V-3 
2JV 


N-l 




N-l 1 
2 l 


2 
V 


N 




N-l 
2 


1 



We may calculate the probability distribution P(E min ; M, N) of the minimum energy for 
M — 1 and 2 by trials of the possible variants. In the case of a single impurity (M = 1), 



21 



there is a probability N+/N = (N + 1)/2N for this impurity to fall on the f-bond, after 
which the minimum energy increases by unity: P(E = —N + + 1) = (N + 1)/2N There 
is also the probability N_/N = (N — 1)/2N for the impurity to fall on the a-bond, after 
which the minimum energy remains unchanged: P(E = — (iV + l)/2) = (N — l)/N Similar 
considerations can be conducted for M = 2. Let C = y/N denote the number of spin clusters 
with the same q on a lattice without impurities. The results of this analysis are summarized 
in Table 2, where each line indicates (I) the number of impurities M, (II) the coupling 
constant of a bond on which the impurity falls (with the signs "+" and "— " corresponding 
to bij = land —1, respectively), (III) the minimum energy E min and (IV) the probability 
P(E min ] M) of obtaining a given minimum energy value. 

With p denoting the probability that a bond contains the impurity J' = —1, the prob- 
ability to find M impurities is P(M;p) = M \^L M y P M (1 — p) N ~ M ■ Thus, a mean value of 
the minimum energy E in the Potts model is 

N /VI 

= J9 M*(N — M)! - P » P ^ M ' N) < A1) 

Calculating the first three terms for this sum using the data in Table 2, we obtain an 
approximate formula for the mean minimum energy 



(E mrn ) ~ (1 - p) N - 2 ((1 - p) 2 P(E; M = 0) + Np(l - p)P(E; M = 1) + 

p 2 P{E; M = 2))+ o{p 2 ) { } 



This formula is applicable if the probability for three impurities to appear on the lattice is 
small, that is, if ^N(N — 1)(N — 2)(1 —p s j N ^ 3 p 3 <^ 1 Upon calculating the probability of the 
impurity concentration for M = 1 and 2 and using expression ( |A1|) , we obtain a formula for 
(E m i n ) at small p. For example, in the case of M = 1 this yields 

/ iV + l\/iV-l\ / JV+1 WiV+l\ JV+1 iV + 1 
(E mm (M = l))= — __) + ( — + 1 —— = — + 



2N ) V 2 ) \ 2N J 2 2N 

By the same token, for M = 2 we obtain 

/w «nx iV + 1 N + 1 12(7-8 
(E min (M = 2)) = — + 



N N(N-l) 



Using these expressions and taking into account that C = \/N , we obtain an expression 

N 



for the mean normalized minimum energy < e min >= (-^f 1 ) as a function of the impurity 
concentration p (in the case of p <C 1) 



< e min (p) > 



j_ 

N 



1 - p) N P(E; M = 0) + Np(l - p) N - 1 P(E; M = 1) + 



-N(N - l)p 2 (l - p) N ' 2 P(n- M = 2) + o(p 2 
1 A 1 \ 1 A 3 \ 26^-4 2 ^ 



(A3) 



2\ 1 -N +P 2 { 1 -Nr P ^^ + 0iP 



22 



Analogous calculations for the concentrations 1 — p <C 1 yield 



(Er, 



N 
N 



1) ) 

2) ) 



JV-1 

'n 2 -i 



N-5 



i$-5 IOC— 14 



iV 



N(N-l) 



and 



p 7V P(E; M = JV) + iV(l - p)p N ~ l P(E; M = N - 1) + 



1 



^N(N - 1)(1 - p) 2 p N ~ 2 P(n; M = N - 2) + o(p 2 ) 



2 V A^ v ^ ; 2 V 



:i-p) 2 ^4r^ + o((i-?) 2 ) 



A' 



(A4) 



The results of calculations using the formulas presented in Table 2 and the results of 
numerical analysis for q = 4, L v = 3, and Lh = 6 (corresponding to a knot with a btimesb 



square lattice diagram) are plotted in Fig.|l2|a in coordinates of the minimum energy modulus 
versus probability. Figure |T2|a shows the results for a fixed number of impurities M — 1 
(squares) and 2 (circles), while Fig.|i~2]b presents data for the fixed impurity con-centrations 
p = 0.005 (crosses), 0.01 (squares), and 0.015 (circles). The results of calculations using the 
formulas from Table 2 and relationship ( |A2| ) are depicted by lines and the values obtained by 
the Monte Carlo method are represented by symbols. The data were averaged over A^ = 10 5 
realizations; the statistical error for most of the points is smaller than the size of symbols. 
As can be seen, the numerical data remarkably fit to the analytical curves. In |l2|b, the 
analytical calculations are illustrated only in the range |-Emm| = H -S- 15, where the results 
could be obtained by expansion into series with terms on the order of p 2 . Thus, we have 
numerically verified the results of analytical calculations performed in this section. 

Let us make some remarks concerning the form of the normalized minimum energy 



< e min (p;N) >- 

e min (p;N) >= 
troughs. As A^ 



N 

in(p\N)\ 



N 



for A^ — > oo. As can be seen in Fig.|9|b, the plots of < 
for the lattices with A^ = 25 and 49 appear as asymmetric 



oo, the values of this function at p = and 1 



< e min (p = 0;N) >= 0.5 1 



and 



<e min (p=l;N) >= 0.5(1 - 



1 

N 
1 

N 

will be equal and the profiles will be symmetric with respect to the transformation p — > 1 —p. 
The bottom will occur on a level of < e(p = 0.5; A^ — > oo) >= 0.415(7) (see Fig.|TT|b). As can 
be seen from formulas ( |A3| ) and ( |A4p , the second derivative of this function with respect 
to p at the points p = and 1 turns zero for A^ — >. Probably, all derivatives of the higher 
orders will turn zero as well. In this case, the derivative exhibits a break at the points 
corre-sponding to the trough "corners." 
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Figure Captions 



Fig.l. a) A knot diagram on the JV = 3 x 3 lattice and b) the diagram splitting. Open circles 
indicate the spin positions in the Potts model; dashed lines show the graphs on the Potts lattice. 

Fig.2. The Reidemeister moves I, II, and II. 

Fig.3. a) A knot on JV = 5 x 5 lattice; b) Potts spin configuration corresponding to this lattice, 
and c) the same spin configuration reduced to a rectangular lattice x L v = 3 x 6. 

Fig.4. A schematic diagram illustrating the step-by-step generation of subdivisions into clusters 
for a column of spins corresponding to the knot state vectors. 

Fig.5. a) An example of the configuration of bonds in which a minimum energy of the Potts 
model is not reached for q = 2; b), c) the arrangement of f- and a-bonds on the Potts lattice with 
Af = 3 x 3 and the impurity concentrations p = and p = 1, respectively; (c, d) the arrangement 
of bonds between col- umns with even and odd numbers, respectively. 



Fig.6. a) A joint probability distribution P ^f £L j of the normalized maximum polynomial 

power jj and the normalized minimum energy E ™ in described by the level curves with a step of 
0.001 for a lattice with A = 49 and an impurity concentration of p = 0.5; b) The coefficient of 
correlation between the maximum power of the polynomial invariant and the minimum energy as 
a function of the lattice size L; dashed curve shows the approximation of these data by a power 
function of L. 

Fig. 7. The distributions of a) probability P(^;p) of the normalized maximum polynomial power 
% and b) probability P (^%^;p) of the normalized minimum energy modulus P m in for a 7 x 7 
lattice and various impurity concentrations p = 0.1, 0.2 and 0.5. 

Fig. 8. The probability of a) the normalized maximum polynomial power and b) the normalized 
minimum energy modulus 1 for the square lattices with N = 9, 25, 49, and 81 and an impurity 
concentration of p = 0.5. 

Fig. 9. Plots of a) the mean normalized maximum polynomial power (-^) and b) the dispersion 
W% not (p) of the value distribution versus the impurity concentration p in the range from to 



3, 5, 7, and 9 and an impurity concentration of p = 0.5; dashed curve in b) shows the results of 
approximation by the power function 0.4185(7) + 0.119(3)L _L11 ( 4 ) . 

Fig. 12. A comparison of the results obtained by analytical (anal.) and numerical (num.) methods 
for the minimum energy probability distributions on a lattice with N = 5 x 5: a) P(\E m i n \; M) 
with M = 1 and 2; b) P(\E min \;p) with p = 0.005,0.01, and 0.015. 
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